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ABSTRACT 


The objective of this study is to assess quantitatively 
the process of "self-development" (Sutcliffe and Forsdyke, 
1950) in extratropical cyclogenesis. The basic approach is 
to determine within the framework of the quaSi-geostrophic 
equations the re eae of thermal advection and latent 
heat release to changes occurring in the intensification 
rate of a surface cyclone which result from modifying the 
field of vorticity advection. The relevant concepts and 
equations are applied to a major cyclone over the central 
United States and to a structurally Simple analytic model of 
meeyeclonic disturbance for which varying horizontal and 
vertical distributions of latent heat releaSe are prescribed. 

Results of computations utilizing the analytic model 
indicate that the effect of latent heat release in self- 
development can be comparable with that of thermal auc con: 
The Separate and, more so, the collective influence of these 
mechanisms upon acceleration of the deepening rate are 
considered synoptically significant; however, the influence 
of translation of the vorticity maximum aloft with respect 
to the surface system appears to dominate changes in the rate 
of development. The quantitative reliability of the results 
obtained for the real data case are in doubt because of 
apparent difficulties in the objective scheme (HOBAN) employed 


ieanalysis of input data. 
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i See T RODUEGT ION 


Pe) OBJECTIVE 

Significant weather in middle-latitudes is almost 
invariably associated with the development of extratropical 
eyclones. ior Surprisingly, therefore, the phySical and 
Gynamical mechanisms responsible for cyclogenesis have been 
subjects of considerable inquiry. The object of such inguiry 
has been to understand more fully the workings of extra- 
tropical cyclones, formulate their governing laws, and 
accurately predict or even control the associated weather. 
Since there are several excellent sources? which summarize 
and synthesize existing concepts of cyclogenesis in middle- 
latitudes, no attempt will be made here to review the subject. 

This thesis concerns the particular aspect of extra~- 
tropical cyclogenesis referred to by Sutcliffe and Forsdyke 
(1950) as "self-development." Specifically, the objective 
of this investigation is to asSessS quantitatively the 
importance of the self-development process and determine the 
relative influence of the contributing mechanisms. Although, 
qualitatively, self-development is readily comprehensible 
and Seemingly important, quantitative assessment of this 
process, to the best knowledge of this investigator, has not 


meen documented in the literature. 


log particular note are the books of Petterssen (1956), 
Palmen and Newton (1969), and the thesis of Downey (1972). 





B. SELEF-DEVELOPMENT 

As discussed by Palmen and Newton (1969), self-development 
is the process wherein progressive deformation of the upper 
current in the vicinity of a cyclone, as affected by changes 
in the thermal structure of the troposphere, results ina 
progresSive increase of the upper-level vorticity advection 
and consequent intensified development of the surface low. 
The process is schematically illustrated in Fig. l. 
Subsequent to the initiation of cyclogenesis, as the cyclone 
achieves a significant circulation, warm advection ahead and 
cold advection behind the developing low center deforms the 
thermal field in the manner shown by Fig. 1. Deformation of 
the mean isotherms (i.e., thickness lines) through the layer 
considered can be identified hydrostatically with the growing 
amplitude of the upper-level contour pattern. The increasing 
amplitude of the upper wave corresponds to an enhanced 
variation of the vorticity (as expressed by the curvature) 
and, therefore, to an increase of the vorticity advection 
along the upper current over the surface low. This, in turn, 
implies an increase in the upper-level divergence, which 
favors progresSively stronger low-level convergence and 
intensification of the cyclone. Note that at each instant 
of time thermal advection primarily contributes directly 
only to motion, rather than to development of the sea-level 
(Or 1000 mb) system. That is, at a given instant, the 
Pattern Of thermal advection acts primarily to create low- 


level cyclonic vorticity ahead of the low and anticyclonic 
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wMemtacity behind it, thereby accelerating the low along a 
path directed from the region of maximum cold to maximum 
warm advection. At that same time, though, the thermal 
advection tends to enhance the primary deepening mechanism,” 
namely the vorticity advection over the surface low. 

If only horizontal temperature advection were important 
in producing thermal changes in the troposphere, an unbounded 
deformation of the thermal field and upper flow would he 
implied. Unbounded increases in amplitude of the upper wave 
does not occur, however, because adiabatic cooling and 
warming associated with vertical motions dampen the influence 
of thermal changes associated with temperature advection. In 
other words, cooling accompanying aScent and warming 
accompanying descent in the regions of warm and cold advection, 
respectively, tend to counteract the influence of horizontal 
thermal advection. Since the magnitude of vertical motions 
increase as the developing circulation intensifies, the 
influence of thermal advection is eventually balanced by the 
opposing influence of vertical motion. Thus, the process of 
self-development calls into play a "Self-limiting" process 
that prevents development from proceeding without bound. 

Mieadeielon to the influence of thermal advection and 
vertical motion, thermal changes in the troposphere can be 


produced by condensational heating. Precipitation and the 


aoaascion here is exclusive of the direct influence of 
latent heat release. 


deel: 





concomitant release of latent heat as, for example, might 
occur ahead of the surface low center and in advance of the 
warm front would result in greater amplification of the 
upper ridge downstream from the low than would be the case 
if warm advection and adiabatic cooling were the only 
processes operative. Self-development would then proceed 
further and the cyclone achieve greater intensity than if 


precipitation were absent. 


eee, BASIC APPROACH 

As noted, the goal of this investigation is to quanti- 
tatively assess the process of dell Fado el oon Thesoacwe 
approach is to determine, within the framework of the gquasi-~ 
geostrophic equations, the magnitude and relative influence 
of the contributing mechanisms. More specifically, 
determination is made of the contributions of thermal 
advection and latent heat releaSe to changes occurring in 
the intensification (i.e., Aeeanamne) rate of a surface 
cyclone which result from modifying the field of vorticity 
advection. > The relevant concepts and computations are 
applied to a structurally simple analytic model of cyclonic 
and anticyclonic disturbances in the baroclinic westerlies 
(after Sanders, 1971) and to a major cyclone over the central 


Wmtted States, 1200 GMT 10 December 1971. 


3The effect of vertical motions within the quasi- 
geostrophic framework appears implicitly in tempering the 
influence of thermal advection and latent heat release, rather 
than explicitly as a separate contributing mechanism. 


| 
NO 





Pee seoRMULATION OF THE PROBLEM 


a a a a te) 


Pee QUASI-GEOCSTROPHIC FRAMEWORK 
It 1S assumed that atmospheric flow is governed by the 


gquaSi-geostrophic vorticity equation 


Q 
SN 


0 W 
cic ae — mee fe + + Pails 
rE ONE ane ty 3p (1) 
and by the thermodynamic equation 
a (8%) = _y.y 2¢ - mee 
ye (ay) VV Xp o(p)w i: H , (2) 


where 7 1S the relative vorticity; w is the vertical velocity, 


dap, 
ati to 


ce temene GeOpotential tendency; Ris the gas consant for dry 


1S a conStant value of the coriolis parameter; 


Slay > C,, is the specific heat for dry air at constant pressure; 


H = ce is the diabatic heating; and V = oK x Vd is the 
O 

geostrophic wind. 

When the geostrophic relationship, % = = V6 , LS used 

O 
to express the vorticity, (1) becomes 
2G INR pene eee ae - 2 dw 

V ra ie V CG aeck = Te (3) 

fombination of (2) and (3) in order to eliminate oe yields 


the quasi-geostrophic omega equation (see, e.g., Holton, 1973), 


le: 











; (1) 


2 2 2 Lc: Kee 
(oV ue ee ©. tee a ee 
(2) 
~ yriyey oo 
Vat vey Tp) 
(3) 
- SS Vn : (4) 
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The solution of this equation yields values of w and, 
hence, = witeneenanpbe SOlution Of (3) for the field of 
geopotential tendency. 

The quaSi~geostrophic equations in the above or related 
forms have been utilized widely (e.g., Sanders, 1971; Stuart, 
moet Danard, 1964) in determination of fields of vertical 
motion and/or geopotential tendency associated with the 
forcing functions on the right-hand side of (4): differential 
vorticity advection, term 1, Laplacian of thermal advection, 
femme 2, and Laplacian of diabatic heating, term 3. In this 
investigation emphasis is not to determine the direct influence 
of the forcing processes upon the fields of vertical motion 
Or surface geopotential, but rather to evaluate their effect 
Upon self-development. That is, attention is focused upon 
the influence of the respective mechanisms upon the field of 
vorticity advection and the resultant changes that would be 
expected in Eee ae eon rate of a surface cyclone. 

Pepliuicarcton OL: (lL) to a leyel near the ground at the 


center of a cyclone (i.e., at the vorticity maximum where 
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V(ttf£) = 0) yields the following expression for surface 
development in terms of the rate of vorticity production by 


mew-ievel convergence: 


oe ~“C Mose: (5) 


Vee ee (6) 


Here, = is the local change in a coordinate system moving 
eee the vorticity maximum. 

Application of the operator, ec BOp CO) yeelds an 
expression for the rate at which the intensification rate 


varies in time aS a conSequence of changes Occurring in the 


Pred Of vertical motion: 


Or 


2.0 CRore ae 6 Ww) 


Temporal variations in the field of vertical motion result 
from changes Occurring in each of the forcing functions of 
(4). Of interest in self-development are variations in the 


weztical motion over the low center and, therefore, in the 





Gbawe On Suritace development that result from changes in the 
meld Of vorticity adhesion. 
In order to evaluate the mechanisms which produce changes 
an the field of vertical motion by modifying the pattern of 
VWOomticity advection, the — Operator is applied to the left- 
hand side and term (1) of (4). With the assumption that the 


stability factor, o, is invariant in time, the resulting 


expression is 


2 (A) 
Bb: d Owes 0 OV. lee 
(oV * £5 5-7) 3 ~ to tp | BG aie) | 
p O 
p 
é IL z 
4 te ap [v-V Se tee a) (8) 


It can be seen from (8) that temporal variations in the 
field of vertical motion result from altering the advecting 
velocity, term A, and the gradient of relative vorticity, 
term B. Modification of each results from changes occurring 
in the field of geopotential, which can be obtained from 


mammeerom Of (4) and Substitution of the results into (3). 


an a developing cyclone the Laplacian of thermal advection 
over the low center is characteristically small in magnitude 
(see e.g., Petterssen, 1956). Time variations in the pattern 
of thermal advéction are therefore not likely to contribute 
Significantly towards modifying the rate of development. On 
the other hand, although precipitation is typically absent 
ma@om the center Of storms in the initiation of cyclogenesis, 
the Laplacian of diabatic heating over the surface low need 
not necessarily be small. Temporal variation in the amount of 
mnrall being produced, for example, by a squall line along 
the cold front can therefore alter the rate of a storm's 
femelopmerce sThis effect, as true aiso for time variation in 
the pattern of thermal advection, is not conceptually included 
in the process of self-development and is not here a subject 
Of investigation. 








The values of 2 Socaimea from (8) on the basis of the 
computed geopotential tendencies can then be utilized to 
evaluate the acceleration of the intensification rate, 
oo (y), from (7). Note that one can obtain the separate 
Semeriobutlons Of vorticity advection, thermal advection, 

and diabatic heating to terms A and B of (8) and, therefore, 
the separate contributions of each to a (y) through the 
effects represented by terms A and B. Actually, the influence 
of term A is not specifically included in the concept of 
self-development, nor 1s the effect of vorticity advection 

on term B. The influence of each, however, is evaluated in 
this study. Figure 2 summarizes the mechanisms contributing 


to time changes of the intensification rate which are evaluated 


icons study. 


B. MODEL AND PROCEDURES 

The concepts just described for accessing the acceleration 
of the intensification rate were applied to a structurally 
simple analytic model of cyclonic and anticyclonic disturbances 
in the baroclinic westerlies and to a major cyclogenesis 
Seer the central United States, 10 December 1971.. In solution 
of the relevant eguations, a three-dimensional grid of data 
points is employed. The vertical grid array (Fig. 3) employs 
Bacal at Gach LLOQ millibar level from 1000 to 100 mb. The 
horizontal grid measures 27 x 27 points with a grid interval 


erel60 km at 60° and 30°N on a Lambert Conic Projection 


Le 





(Seeunigss3 £or grid location). All derivatives are 

replaced by central finite differences, with the resulting 
set of algebraic equations solved by the Liebmann relaxation 
process. The boundary conditions are w = 0 and oe ae (tcc Glilah 
lateral and vertical boundaries. The static stability 
factors utilized are those for the standard atmosphere as 
fieem in Table ~, The standard value of the coriolis factor, 
fees) that for 39°N, .92 x Weecece % 

As shown in Fig. 5, subSequent to generating the three- 
dimensional fields of geopotential and heating, the first 
step in the computation procedure is to evaluate the contri- 
bution of the forcing functions in (4) to the field of vertical 
motion, w. The respective contribution of each to the 
geopotential tendency is then computed via (3). Actually, two 
sets of tendencies for each forcing function are obtained; 
the first being that associated solely with the divergence 
(38) term of (3), while the second is the total tendency. The 
first set of tendencies relates to the developmental influence 
of each forcing function, while the second set includes the 
effects of development and motion. 

Once the field of geopotential tendency is obtained, the 
influence of each contributing factor thereto upon = 1s 
evaluated from (8). The separate contributions of terms A 
and B are determined. Finally, the respective influence of 


the relevant physical effects to temporal variations in the 


rate of development is ascertained from (7). 
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The forcing functions of (4) and , consequently, all 
which follows therefrom are specified by the three-dimensional 
field of geopotential and distribution of diabatic heating. 
Input for the real case study was height data at mandatory 
Gata levels (1000,850,700,500,400,300,200, and 100 mb), 
converted by quadratic fit where necessary to give ¢ at each 
100 mb level. Analysis of the geopotential fields were 
accomplished by an objective technique known as HOBAN 
(Inman, 1970). HOBAN is basically a successive approximation 
Scheme that starts with a large influence radii and reduces 
the radii on succesSive passes to give a more local fit on 
the last pass. The scheme makes use of height and wind 
within and outSide a local limited area grid. The scheme 
employs a distance dependent weighting function following 
Cressman (1959). In an extensive evaluation of HOBAN, 

Stuart (1974) found the analyzed height field and quasi- 
geostrophic omegas derived therefrom compared quite favorably 
with the height and omega fields obtained from carefully 
prepared and consistent subjective analyses that were used 

aS the basis for judgement. Analyses of the 1000 and 500 mb 
height fields on 1200 GMT 10 December 1971, the mane fOr 
which computations are performed, are shown in Fig. 6 and 

7, respectively. 

In addition to the real data case, computations are 
performed for the structurally simple model of the geopotential 


and temperature fields of cyclonic and anticyclonic 











aeetimirbances, aLter Sanders (1971). At 1000 mb, the 


disturbances of geopotential iS given analytically by 


cau 
L 


27 


Cos (Gc) Cos Sree (9) 


d(x,y,1000) = 1000 


where the x and y axes are directed eastward and northward, 
respectively; byo00 is the amplitude of the pressure 
perturbation; L is the horizontal wavelength; and A is the 
phase lag of the 1000 mb geopotential field relative to the 
temperature field, which is specified by 

fis, yY,p) = Ty (P) -~ (l-aln =S2") (ay + TCos of x. CoS £m y) 


(10) 


Here, Ty (P) is the mean temperature at the pressure level, 

p, “a" is the mean north-south temperature gradient, and 

. is the amplitude of the temperature perturbation. The 
vertical structure of the temperature perturbation is governed 
by the factor (1 - aln =) which allows for damping with 
elevation and reversal of sign at a pressure level determined 
by the value of a. With a = .722, the temperature reversal 
and, hence, the simulated tropopause occur at the 250 mb 
level, as in middle-latitudes in the real atmosphere, while 
the strength of the temperature perturbation at 506 mb drops 
memone-nalt its 1000 mb value. Note that A = 0 corresponds 


to a system of warm lows and warm highs. The 1000 mb pattern 


with A = L/2 yields cold lows and warm highs. The 1000 mb 


a0) 





pattern with A = L/4, representing the typical intensifying 
Situation in the real atmosphere, and with model parameters 
obtained from the system of 1200 GMT 10 December 1971 is 
shown in Fig. 8. 

At all levels above 1000 mb, the geopotential fields are 


obtained by hydrostatic integration; thus, 


p 
¢(x,y,p) = 6(x,y,1000) + sf $8 (x,y,p) ap 
Looe 3 
oo8 7 1000 _ a 1000, 2 
= by (P) + Oe ayn Oe) [R 10a “> (R 5) (1n et ] 
me (avetel Cos ue Cos ony jaa; (11) 
where 
1000 
d4(P) = R i T,(p) din p eles 


The height field at 500 mb corresponding to the Situation in 

Fig. 8 is presented in Fig. 9. It should be noted that 

Fig. 8 and Fig. 9 represent contour analyses of the grid 

point height values ascertained from the analytic expressions 

above. Further, all computations are performed on the 

twee 27 x 10 grid in the same fashion as for the real caSe. 
The main contribution to diabatic heating 1s assumed to 

be from the heat released by condensation. Other heat sources 

and sinks are neglected. Furthermore, it is assumed that 

condensation occurs as liquid water which immediately falls 


Out as precipitation. 


AL 





For the real data case, precipitation rates at 1200 GMT 
10 December 1971 were obtained from the observed two-hour 
amounts ending at 1300 GMT, subjectively averaged over 
160 km? areas centered at grid points (Fig. 10). Area 
averaged precipitation rates were converted to the vertically 
integrated heating by multiplication with the latent heat 
factor, L. The vertical partitioning of the latent heat 
release is per the analytic model discussed below, with the 
level of maximum heating, 400 mb, consistent with the fact 
that the observed precipitation was primarily convective in 
nature. 

In the analytically specified cyclonic system, the baSic 
horizontal distribution of rainfall is modeled as an ellipse 
in which the maximum precipitation rate, Pays 1s at the center 
and values of P decrease exponentially therefrom (after 


iereton, 1973). ‘Thus 


NO NA), 7 Py e f 


where A and B are scale factors which specify the major and 
minor axis (x and y, respectively) of the ellipse defined 
pee x,y) = 1 Pia The vertical distribution of heating is 


Srecne form 


e —*4 sin z , (1-3) 
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where 


1000 mb 


1G: 

! 

Ee. 
— 

16 

[I 


200 mb = level at which heating 
GoOesTcO ZErO 


O 
bem 
II 


and c specifies the level of maximum heating. 

Computation of all quantities are performed for four 
Melations Of the basic distribution of heating. The first, 
Case I, is that of an ellipse centered one grid to the east 
and One grid to the north of the 1000 mb low, with the value 
of A and B such that the major and minor axes are 1440 and 
640 km respectively (Fig. 11). The value of “c" in this 
case 1S in accordance with the heating maximum in the vertical 
ieeeag at 7/00 mb (Fig. 12). Motivation for the first case is 
to Simulate the stratiform precipitation characteristically 
observed ahead of a low in advance of the warm front. The 


maximum intensity of heating at 700 mb is explained by large 
dq 
dp 
Surface and large values of vertical motion near 500 mb. 


dq 
The correlation between w and oa , which determines the 


rate of stable precipitation, therefore tends to be largest 


S 





values of (qo = Saturation specific humidity) near the 


somewhere between these levels. 

The second case considered, Case II, is the same as the 
first, except that the level of maximum heating is specified 
as 400 mb (Figs. 11 and 12). Here, an attempt is made to 


Simulate the effect of convection embedded in a warm-frontal 
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precipitation pattern. As noted by Tracton (1968), warm- 
front precipitation is largely produced by convective 
S erturning, rather than by stable saturated ascent. Although 
little is known about the detailed distribution of latent 
heat release by cumulus convection, theoretical treatments 
Sreevarious investigations (e.g., Kuo, 1965; Kasahara and 
Asai, 1967) appear to agree that convective heating is 
largest in the middle to upper troposphere; hence, the choice 
@eme00 mdb for Case II. 

ine third variation of heating for which computations 
are performed, Case III, is that of an elongated ellipse, 
1440 x 320 km, centered one grid distance to the east and 
one grid distance to the south of the 1000 mb low and oriented 
Seeetrom the north-south direction (Fig. 13). Motivation 
in this case is to Simulate the influence of asquall line 
ememg Or jJuSt ahead of the cold front. Squall~line precip- 
jMeseron 15, Of course, convectively produced and, hence, 
the level of maximum heating is specified to be at 400 mb. 

The final case examined, Case IV, represents a super- 
position of the squall line (Case III) and the stratiform, 
warm frontal precipitation pattern (Case I). 

Note that in all cases, P,, 1s specified to have the 


M 


synoptically reasonable value of .5 in hr + 
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Pie scou lho Debi oCusSs LON 


PeeeONALYTIC MODEL CASES 

Table II presents the tendencies at the 900 mb low 
peeer.- The tendencies are, of course, a measure of the 
instantaneous rate of low-level development but are presented 
in synoptically familiar terms of change in central pressure 
per hour. Note that the tendencies associated with 
vorticity advection, Xe and thermal advection, Xo, are 
independent of the diabatic heating, X31 and, therefore, 
the same for each of the cases considered. 

Exclusive of the influence of latent heat release, it is 
clear from Table II that vorticity advection is the dominant 
process contributing to development, i.e., the active 
deepening mechanism. Thermal advection, in accord with the 
Qualitative arguments of Petterssen (1956), produces pressure 
falls in advance and rises to the rear of a cyclone, but has 
very little effect at the low center itself (Fig. 15). The 
deepening produced by vorticity advection, -.63 mb hr -, is 
consistent with the value obtained by Sanders (1971). 

The influence of latent heat release upon the deepening 
rate varies according to the distribution of the heating with 


respect to the low center. Values range from -.55 mb hr + in 


= 

~The 900 mb level is the lowest one for which tendencies 
are explicitly computed. Tendencies at 1000 mb required for 
boundary values in solution of (8) are obtained by quadratic 
Mire rpOolation. 


a5 


| 





the squall line Simulation (Case III) to -2.36 mb hr iByg 


the case where the level of maximum heating associated with 
warm frontal precipitation is at 700 mb (Case I). Of course, 
in each of the cases considered a shift in the location of 
the center of heaviest precipitation would affect the 
magnitude of the computed influence of heating upon the 
Beepening rate. Translation of the maximum in precipitaticn 
[mento the low center, aS 1S characteristic of the initial 
stages of cyclogenesis (Tracton, 1973), would result in 
increasing the influence of condensational heating upon the 
instantaneous rate of development (e.g., See Fig. 16). Note 
also that the magnitude of the tendencies associated with 
heating are proportional to the specified values of Pare which 
in all the above cases is .5 in keane 

Values of the computed time rate of change in intensifica- 
tion rates at the 900 mb low center appear in Tables III and 
Iv. The results presented in Table III are based upon 
geopotential tendencies associated with only the divergence 
term of (3), while Table IV reflects the field of total 
tendency produced by each forcing function. Reference should 
be made to Fig. 2 for the meaning of terms in Tables Ti send 
Iv. Note that calculated values of the instantaneous 
acceleration of the intensification rate have been applied 
to 12 hour periods for ease in interpreting the results. 


The terms which specifically relate to the process of 


OX OX 
2B 3B 
re and ae Tie fLorier Certreccs 





self-development are 
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the influence of thermal advection upon the gradient of 
relative vorticity and, hence, vorticity advection over the 
storm center, while the latter represents the effect upon 
the gradient of relative. vorticity resulting from diabatic 
heating. Table III indicates that the thermal advection 
effect produces an increase in the deepening rate of .46 mb 


hr 


per 12 hr. This corresponds to an increase of 75 per 
cent in the rate of development associated with the field 

of vorticity advection. The physical mechanisms which are 
operative here are illustrated in Figs. 17 to 19. In 
response to the pattern of thermal advection (Fig. 17), height 
rises occur aloft in advance of the surface low, while height 
falls occur aloft behind the low position (Fig. 18). The 
height rises and falis correspond to generation of anti- 
cyclonic and cyclonic vorticity, respectively. The enhanced 
Variation of vorticity and resultant increase in vorticity 
advection along the upper current over the surface low 
intensifies the vertical motion over the low (Fig. 19). The 
consequent increase in low-level convergence then results in 
an increased rate of development. 

The influence of latent heat release in self-development, 
as shown in Table III, ranges from -.002 to -.50 mb hr per 
12 hr. The smallest effect occurs in the squall line 
Simulation (Case III). This reflects the fact that the 
height changes induced aloft by the heating result in an 


increased gradient of vorticity across, rather than along 


the upper flow. Hence, the field of vorticity advection 
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over the low is not Significantly altered. In addition, 
as can be seen in Fig. 19, the impact of the heating upon 
the vertical motion 1S principally at levels above 700 mb. 
Consequently the effect of heating upon the low-level 
convergence and, therefore, upon self-development is minimal. 

The largest impact of latent heat release occurs in the 
700 mb warm front precipitation and in the combination 700 
mb warm frontal and squall-line precipitation case (Cases I 
and IV, respectively). Note that the computed acceleration 
of the intensification rate, -.50 mb hr eyeic LZ Ages “shins eyenayel 
cases, 1S comparable in value to that ascertained for thermal 
advection. Here, the upper level height rises induced by the 
heating (Figs. 20 and 21) tend to amplify the ridge downstream 
from the low center, such that the enhanced vorticity 
advection over the low reSultsS in an increased rate of 
intensification. The lesser effect (-.035 mb hie per 12 hr) 
when the level of maximum heating is at 400 mb, rather than 
700 mb, results from the shift upward of the level at which 
heating has its greatest influence on the vertical motion. 
The impact of the heating in Case IV is Simply not felt at 
low levels. | 

Recall that the motivation for the 400 mb warm frontal 
precipitation case was to simulate the effect of convection 
embedded in large-scale stratiformprecipitation, which alone 
would have a level of maximum heating around 700 mb. Actually, 


in most inStances, Said convection does not extend above the 
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ago to 500 mb level (Tracton, 1968), so that the effective 


level of maximum heating would likely lie Somewhere between 


OXa55 


ae 
frontal precipitation is thus -.50 to -.035 mb hr - per 12 hr. 


700 and 400 mb. The total range of values of for warm 
While the apparent influence of convection is to reduce the 
effect of heating upon self-development through raising the 
level of maximum heating, it should be noted that the 
maximum in the rainfall rate (to which all computed quantities 
are proportional) is likely greater in the case of convection 
embedded in Stratiform precipitation than when the precip- 
meateron iS purely stratiform. in addition, the horizontal 
scale of heating is less in convective Situations, so that 
through the influence of the Laplacian in the relevant 
equations, the importance of heating would be greater. 

With reference again to Case IV, note that, although the 
combined maximum precipitation rate 1S somewhat greater than 
that of Case I (Fig. 14 and 11, respectively), the effect of 
heating on self-development is exactly the same. This 
reflects the offsetting influence of an effective rise in 
the level of maximum heating and the shift in the location 
@ieeene center of precipitation with respect to the ou 
Selcer . 

eon table £11 16 can be seen that the effect of 
Sereicity advection upon the gradient of vorticity and, 


Guemerore, upon the field of vorticity advection itself 
0X 
1B i 


— oleien dl emabe albs' 


), is to produce an increase of .22 mb hr. 


ZS 





the rate of development. As stated previously, this 
mechanism is not conceptually included in self-development. 
The effects of thermal advection, vorticity advection, and 
feaeing upon the pattern of vorticity advection through 
their influence upon the advecting velocity (all "A" terms) 
can be seen from Table III to be negligible. These effects 
are also not specifically encompassed in self~development. 
It can be seen from Table IV that the ants computed on 
the basis of the total tendency produced by each forcing 


function are either negligible in value or positive in Sign. 
IXop 
7 


Preauces a decrease in the deepening rate of .73 mb hr 





Thus, for example, the thermal advection effect, 
a: 


per 12 hr, rather than the previously discussed increase of 


mao mb hr 


per 12 hr. The apparent reason for this is that 
the structure of the analytical model is such that the maximum 
in the field of vorticity advection lies ahead of the surface 
low position (Fig. 22). Consequently, when one includes the 
vorticity advection term of (3) into the computed tendencies, 
the induced translation of the upper trough eastward results 
in smaller values of vorticity advection becoming superimposed 
over the surface ee thes Lesulmaltes Leauced rate of low- 


level development (.73-(-.46) = 1.19 mb hr ~ pore 2 hr) 


is then sufficien to offset the oppoSing influence of the 


Srhe speed of translation of the upper trough is greater 
jman that Of the Surface low (Sanders, 1971). 
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the developmental components of the forcing functions. It 
should be noted that the Situation here iS similar to that 

of a cyclone in the occluding (but not fully occluded) 

Srage OL its life history. In an earlier stage of development, 
when the maximum of vorticity aloft lies to the rear of the 
storm center, the effect of translation would be to progres- 
Sively increase the rate of development. If the influence 

of translation in Such a case were of the Same magnitude as 
here, this process would dominate the individual contributions 


to self-development. 


Epeeenot OF 1200 GMT 10 DECEMBER 1971 

The synoptic case is that of a developing frontal cyclone 
over the central United States on 10 December i971 (Fig. 6). 
The initial development of the storm occurred around 0000 GMT 
10 December over southeastern Oklahoma in association with an 
outbreak of convective showers in the vicinity of the low 
Senter (Fig. 23). By 1200 GMT, the convection began to 
spread eastward away from the center of the storm as the low 
continued to intensify and move northeastward toward the 
Great Lakes. This association between the initisetion of 
cyclogenesis and convective activity, and the subsequent 
dissociation of the convection from the low center has been 
meted by Tracton (1973) to be characteristic of approximately 
one-half the cases of major cyclogenesis over the eastern 


two-thirds of the United States. 


ow 











Figure 24 presents a plot of the observed change in 
central pressure versus time for the period 0000 GMT 10 
December to 0000 GMT 11 December 1971. The 12-hour average 
Geepening rate at 1200 GMT 10 December, the time for which 
1 


computations are performed, is 1.08 mb hr The increase 


in the rate of development which 1s occurring at this time 


feieG? mb hr — 


Den sh2o hr. 

Before presenting the results for the real data case, 
a few comments are in order with regard to their reliability. 
Although the contour analysis produced by the HOBAN scheme 
at each level appear synoptically acceptable (as, for example, 
Fig. 7) and are generally compatible with the corresponding 
NMC analyses, the vertical consistency between levels is 
somewhat dubious. This is apparent from the irregularities 
in the derived 1000-500 mb thickness pattern (Fig. 25). 
Adjustments to the HOBAN scheme, which are beyond the scope 
of this investigation, will be required in order to resolve 
mmemeroblem. in addition, Stuart (1974) noted that HOBAN 
Goes well on the general pattern and magnitude of main 
Sees, Dut not so well on their precise location. In the 
Gase under study here, this is manifest in placement of the 
center of maximum vorticity advection ahead of the low 
rather than behind it, as is true in the NMC analysis and is 
implied by the Bhset es increase in intensification rate. 

With the above reservations, Tables II, III and IV 
present the computed tendencies and accelerations of the 


tendency at the 900 mb low center. The vorticity advection, 
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Table II, accounts for deepening at the rate of -~.66 mb hr+, 


pies 1S largely offset by the thermal advection term 
whose magnitude is .56 mb he HOBAN analyses do, in fact, 
show cold advection at low levels over the surface low. This, 
in itself, is not surprising in view of the asymmetrical 
nature of the system in question; however, there is little 

or no temperature advection (e.g., at 850 mb) in the 
corresponding NMC analysis. Table II indicates that diabatic 
heating contributes negligibly to the deepening rate, for 

the maximum of the precipitation pattern is too far removed 
from the low center to contribute directly to its development. 


1 


Minestotal computed tendency (-.12 mb hr ~) is clearly 


inadequate to account for the magnitude of development which 
1s observed. 
From Table III one can see that thermal advection 


Semeributes —-.43 mb hoe per 12 hr to the deepening rate 

OX 
Os 
6t 

Givergent effect only). Heating contributes negligibly to 


OXan 
oc 


specifically incorporated in self-development, produces an 


a 


mememgh 1tSs influence upon the gradient of vorticity ( 





Same, ( ), while the vorticity advection, which is not 
acceleration in the tendency of -.43 mb hrs eleva Mie Mone 

The combined effect of the thermal and vorticity advection 
upon changes in the deepening rate account for about one-half 
the observed acceleration; however, one can see from Table III 
that the remaining non self-development terms considerably 


emeset their influence. 
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Metre Grom taole lV that inelusion of the vorticity 
advection term into the tendencies required for computing 
acceleration of the intensification rate renders all computed 
values poSitive. That is, all effects contribute to 
decreasing the rate of development. This stems from the 
erroneous placement by the HOBAN scheme of the maximum of 
vorticity advection ahead, rather than behind the low 


center. 
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ive OUMMARY AND CONCLUSIONS 


The Objective of this investigation was to assess 
quantitatively the process of self-development. Computations 
of relevant quantities within the framework of the guasi- 
geostrophic equations were performed on a real data case and 
an analytic specification of a cyclonic disturbance for which 
varying horizontal and vertical @istributions of latent heat 
release were prescribed. 

The quantitative reliability of the results obtained 
for the real data case are dubious because of apparent 
problems in the objective analysis scheme, HOBAN. Clearly, 
future effort should be directed towards remedying, if 
possible, the difficulties in HOBAN, so that the concepts 
and procedures described herein might be more reliably 
applied to real data case Studies. 

Results of computations utilizing the analytic model 
indicate that when only the divergence (1.e., development) 
term of (3) is utilized in determination of the fields of 
geopotential tendency required for determination of the 
OX 


5¢ 5) the effect of latent heat release in self-development 


IXo5, 


( dt 


) can be comparable with that of thermal advection. 


Specifically, in both the 700 mb warm frontal precipitation 
case (Case I) and the combined 700 mb warm frontal and squall 


line precipitation case (Case IV), the computed time rate 


oo 





of change in deepening rate due to condensational heating 


was -.50 mb hr * 1 


per i2 hr, versus -.43 mb hr ~ per 12 hr 

due to thermal advection. In light of characteristically 
observed accelerations of intensification rate, as, for 
example, the —-1.66 mb hr per IZ hr for the 10 December 1971 
case, the contribution of each mechanisms is considered 
synoptically significant, and their simultaneous occurrence 
even more so. 

When the vorticity advection term of (3) is included in 
determination of the fields of geopotential tendency induced 
by the respective forcing mechanisms, the resulting SA's 
indicate a decreaSe, rather than an increase in the rate of 
development. This stems from the fact that the maximum of 
vorticity advection in the model lies ahead, rather than 
behind, the 1000-mb low center, which iS contrary to the 
Situation in the real atmosphere in aSsociation with a 
progressively intensifying cyclone. The magnitude of this 
effect, which represents translation of the vorticity maximum 
aloft with respect to the surface system, 1S 1.19 mb hr + 
per 12 hr. Clearly, if the magnitude of this influence were 
the same for the Situation where the maximum of vorticity is 
approaching a surface cyclone, it would be the dominant 
meecnanism contributing to intensification of the rate of 
development. 


Note that the influence of heating 1s dependent upon the 


intensity and three-dimensional distribution of precipitation. 
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Mimtenis study, the intensity and distribution are synoptically 
and physically reasonable but, nevertheless, somewhat 
arbitrary. The same is true for all the structural parameters 
of the model, such as wavelength, intensity of the temperature 
perturbation, etc. Further experimentation with the model 

is therefore in order to generalize the results over a wide 
spectrum of situations. Finally, the concepts and approach 
employed here can be an important tool in diagnosing the 
performance of numerical prediction models in predicting 


extratropical cyclogenesis. 
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Table II. Geopotential tendency at 900 mb low center. 
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Table III. Time rate of change of tendency derived from 


divergence term of (3). Units: mb hr 


(- = increase in deepening rate) 
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Table IV. Time rate of change of tendency derived from 
divergence and vorticity advection terms of (3). 


Units: mb hrt per 12 hr (- = increase of 
Geepening rate) 
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Figure 1. Schematic 500-mb contours (heavy solid lines), 
1000-mb contours (thin lines), and 1000-500 mb 
thickness (dashed), illustrating the 
"self-development" process during growth of a 
cyclone. Thermal advection (geostrophic) of 
the 1000-500 mb layer is inversely proportional 
to the size of quadrilaterals formed by the 1000 mb 
contours and 1000-500 mb thickness lines (after 
Palmen and Newton, 1969). 
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Figure 2. Summary of contributing factors to time changes 
of the intensification rate; note that Xon ana X25 
are the only effects specifically encompassed in 


self-development. 
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Figure 3. Vertical grid, prime guantities are obtained 
by quadratic interpolation from data at mandatory 
levels in case of Z' and @6', and from levels 
where guantities are explicitly evaluated in 
case of X', subscript "m" refers to analytic 
specification of initial fields. 


43 











oa 


& 


Ha 
as 4 
Ae 


ho 
~J 
.Q 
ry 
!- 
As 


mequre 4. Placement of horizontal 27 3 


44 





head age CONG vance, erent 
Pe ee ee eee ee won oh ema 
OMe wei: Taomor, COrrol iS, and 


veonotential; enter precipitation 
and distribute vertically 


aft . 2 


vonpute at ae hand side (4) 
Term 1 f Sl v-vlp vss f )] 
Term Z t-V2 (Y- 928 | 
eit ve pV 


Solve (44) for w,, wand Ww, 




























solve (3) aay ance 
divergence alore and total tendency 


eo 


Sompute rish: head side (3) 
Term A 5 [ gy. Vv ( # LV +6)] 














oes 2 
Term B foo [v-v ds “a ae 






Solve (8 


) con Si Sebene $e 
Parr an 


é ot 


rs as Ye a mat oli 












Sole ( ia 


Cian, Pon, Oi 









+ Oia, easier §Xon 


a CR ee 


Figure 5. Schematic illustration of computation procedure 








Figure 6. HOBAN generated 1000 mb height field for 
1200 GMT 10 December 1971. (NMC frontal 
analysis superimposed). Units: m. 
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Figure 7. HOBAN generated 500 mb height field for 
1200 GMT 10 December 1971 (+ indicates 
1000 mb low position). Units: dm. 
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Figure 10. Observed rate of precipitation 1200 GMT 
10 December 1971. Units: in hr?-, 
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Figure 12. Vertical distribution of heating; solid line 


represents heating maximum centered at 400 mb 
and dashed line represents heating maximum 
at 700 mb. 
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Paomre 19. 


©100 


Vertical profile of time rate of change of 
vertical motion, ae , over low center produced 
by thermal advection (solid line) and by 
Giabatic heating of Case III (dashed line). 
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Figure 23. Composite sSurface-radar Summary chart for 
0545 GMT 10 December 1971. 
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Figure 24. Plot of Central Pressure (Units: mb) versus 
time (Units: hrs) from 0000 GMT 10 December 
to 0000 GMT 11 December, 1.971 
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Figure 25. 1000-500 mb thickness pattern for 1200 GMT 
10 December 1971. Units: dm 
(+ indicates 1000 mb low center) 
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